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ABSTRACT 

The detection of Lya emission in star-forming galaxies in different shapes and in- 
tensities (always smaller than predicted for case B recombination) has puzzled the 
astronomical community for more than a decade. Here we use two dimensional cal- 
culations to follow the evolution of superbubbles and of the H n regions generated by 
the output of UV photons from massive stars. We show the impact caused by massive 
star formation in the ISM of different galaxies and we look at the conditions required 
to detect Lya emission from a nuclear H n region, and the variety of profiles that may 
be expected as a function of time. 

Key words: Hn regions - galaxies: starburst -galaxies: superbubbles - ISM: struc- 
ture - ISM: evolution - methods: numerical. 



1 INTRODUCTION 

It has been conjectured that primeval galaxies at large red- 
shift would be easily detected from their Lya emission (Par- 
tridge and Peebles 1967; Meier 1976). Ultraviolet observa- 
tions of nearby starburst galaxies however, have revealed a 
much weaker Lya emission than predicted by simple mod- 
els of galaxy formation (Meier & Terlevich 1981, Hartmann 
et aJ. 1984, 1988, Deharveng et al. 1986, Terlevich et al. 1993, 
Kunth et al. 1997, 1998, 1999). Star-forming galaxies have 
also been recognized at large redshifts and many of them 
show weak Lya emission or none at all (Steidel et al. 1996, 
Lowenthal et al. 1997). The reason for the weakness of the 
Lya emission has been the subject of a series of debates. It 
was early realized that pure extinction by dust would be un- 
able to explain the low observed Lya/H/3 although this was 
disputed by Calzetti & Kinney (1992) who tentatively pro- 
posed that proper extinction laws would correctly match the 
predicted recombination value. Valls-Gabaud (1993) on the 
other hand suggested that ageing starbursts have reduced 
Lya equivalent widths affected by strong underlying stellar 
atmospheric absorptions. Early IUE data has provided evi- 
dence for a correlation of Lya emission and the metallicity 
of the gas as expected if Lya is destroyed by dust absorption 
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(Meier & Terlevich, 1981, Hartmann et al. 1988, Terlevich 
et al. 1993). Dust absorption can be particularly effective if 
multiple scattering in neutral hydrogen selectively increases 
the path length of these photons through the dusty regions 
(Auer, 1968). Moreover high density Hi regions can also 
backscatter photons into the Hll region in which even a 
small amount of dust mixed within the ionized gas could 
be an important agent of destruction. Detailed calculations 
performed by Chen & Neufeld (1994) show how these effects 
can lead to Lya absorptions even for galaxies in which the 
unattenuated spectrum would show a strong Lya emission 
line. 

Then, as also shown by Chariot & Fall (1993) the struc- 
ture of the interstellar medium (porosity and multi-phase 
structure) is most probably an important factor. As a result 
the Lya emission vs metallicity correlation should show a 
large scatter, as evidenced in the recent compilation of 21 
local low metallicity starburst galaxy spectra from the IUE 
archives (Giavalisco et al. 1996). 

New HST data on ten star-forming galaxies (Lequeux 
et al. 1995; Kunth et al. 1997, 1998; Thuan & Izotov 1997) 
has allowed the further realization that the velocity struc- 
ture in the interstellar medium plays also a key role in the 
transfer and escape of Lya photons (Kunth et al. 1997). 
Three types of observed line profiles have been identified: 
pure Lya emission; broad damped Lya absorption centered 
at the wavelength corresponding to the redshift of the Hll 
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emitting gas; and Lya emission with blue shifted absorption 
features, leading in some cases to P Cygni profiles. As noted 
by Kunth et al. (1998), Lya emission with deep absorption 
troughs at the blue side of the Lya profiles, evidence of a 
wide velocity field, has been detected in four of the observed 
galaxies. Interstellar absorption lines (Ol, Sill) are also sig- 
nificantly blueshifted with respect to the Hll gas. 

The velocity and density structure of the neutral gas 
along the line of sight rather than the abundance of dust 
particles alone, seem to be the determining factor for the es- 
cape of Lya photons in these objects. Lequeux et al. (1995) 
in the analysis of the HST data of Haro 2, proposed that 
if most of the neutral gas is outflowing from the ionized re- 
gion, perhaps in a galactic wind resulting from the intense 
star formation activity, then Lya photons will escape (par- 
tially) unaffected. If the Hi is static with respect to Hll, the 
destroyed Lya photons are those emitted by the H II region; 
otherwise, they would correspond to the stellar continuum at 
wavelengths close to Lya. Nine out of the ten galaxies show 
broad Lya absorption, with a very narrow derived range of 
logarithmic column densities (19.7 to 21.5 cm -2 ). 

Being aware of the complications intrinsic to the trans- 
port of Lya lines, we postulate here a simple scenario that 
accounts for the variety of Lya line detections in starburst 
galaxies. The scenario is supported by two dimensional cal- 
culations of multi-supernova remnants evolving in gas-rich 
disk galaxies covering the range between dwarf and massive 
galaxies such as the Milky Way (see Sections 2 and 3). The 
calculations account for the photoionization produced by a 
massive starburst, and its impact on the surrounding gas 
as a function of time. Our models are based on the syn- 
thetic properties of starbursts, as derived by Mas-Hesse & 
Kunth (1991) and Leitherer & Heckman (1995); and con- 
sider a mechanical energy injection rate (Eq ~ 10 38 10 42 

erg s _1 ) which over several 10 7 yr is equivalent to some tens 
of thousand of supernova explosions, and the correspond- 
ing input of ultraviolet photons generated by the massive 
stars in the considered cluster. The continuous mechanical 
energy deposition generates a remnant with an outer shock 
that sweeps the surrounding gas into a thin shell, while the 
evacuated cavity is filled with the hot ejecta; efficiently ther- 
malized at the inner shock. The remnant at first decelerates 
rather rapidly, although its overall expansion allows it to 
depart from spherical symmetry, given the lower densities 
away from the galaxy plane. The latter eventually becomes 
insufficient to cause any further deceleration and the rem- 
nant blows out (see Tenorio-Tagle & Bodenheimer 1988; 
and Silich & Tenorio-Tagle 1998, hereafter referred to as 
Paper I) allowing for the venting of the hot matter into the 
galaxy halo. As shown in Paper I, blowout occurs quite early 
in the evolution, as the superbubbles reach dimensions ~ 1 
kpc. This leads to shell fragmentation and to the venting 
of the processed matter into the extended low density halo 
where the outer shock forms a new shell of swept up mat- 
ter. This in the case of dwarf galaxies, initially expands with 
speeds that well exceed the local escape speed of the galaxy. 
However, its motion into the gaseous halo causes a contin- 
uous deceleration lowering its velocity to values well below 
the escape speed of the galaxy. These facts are in excellent 
agreement with detailed observations of dwarf and irregu- 
lar galaxies which show giant complete remnants expanding 
at a pace comparable or slower that the escape velocity of 



their parent galaxy (see Marlowe et al. 1995, Martin 1998). 
Here we consider the fact that blowout must also allow for 
the escape of a fraction of the UV radiation produced by 
the massive members of the starburst, causing an extended 
H II region into the galaxy halo. The latter is here also held 
responsible for the escape of UV and Lya photons from the 
galaxy. We follow the time evolution of the photoionized 
region (Section 3) accounting for the drop in stellar lumi- 
nosity caused by the ageing of the cluster, as well as for the 
range of recombination times that arises from the different 
densities that characterize the galaxy haloes and the decel- 
erating shells. The calculations allow us to predict (Section 
4) a variety of Lya profiles that are to be expected from 
these sources, as a function of time and inclination angle; 
and a comparison with the existing observations is included 
in the Discussion section. 

2 THE EVOLUTION OF SUPERBUBBLES 
IN GAS-RICH DISK GALAXIES 

2.1 The ISM density distribution 

We approximate our galaxies with the prescription for the 
stellar, dark matter, and gaseous components as given in 
Paper I and in Li & Ikeuchi (1992). The gas density dis- 
tribution allows for two isothermal components, related to 
the central dense molecular core (with a temperature of 100 
K°), and a low density neutral gas halo (with temperature 
1000 K°). 

Pg = pcore+ < Phalo > (1) 

Both components are supported in a quasi-equilibrium state 
by rotation and random gas motions with an effective gas 
pressure 

Pext — — (PcoreCcore-h < Phalo > C^aZo); (2) 

where C CO re is the velocity dispersion of the compact dense 
core, and Chaio that of the extended halo component. 

Typical initial density distributions, for several of the 
cases here considered (see Table 1), are shown in Fig. [l^,- 
c. In all cases the ISM mass amounts to 10 percent of the 
total mass of the galaxy and ranges from 5 x 10 7 M© in the 
case of the smallest dwarf galaxies, to 10 10 M s for the more 
massive disk-like galaxies considered here. 

The initial column densities vary between 10 20 and 10 23 
cm~ 2 for different models. Note however, that only a minor 
part of this value comes from the extended low density halo. 
Most of the undisturbed gas column density is (as shown in 
Fig. ^) associated with the central molecular core. 

2.2 The starburst mechanical energy and 
UV radiation input rates 

The mechanical luminosity from massive starbursts (M s tars 
> 10 6 Mq ) is known to lead to a rapidly evolving superbub- 
ble able to blowout from the gaseous galactic disk configu- 
ration into the extended H I halo, as Eq exceeds the energy 
input rate required for the remnant to reach a typical disk 
scale-height with a supersonic velocity. As shown by Koo 
& McKee (1992) if Eq exceeds the threshold luminosity Lb 
= 7.2 x 10 36 P4-ff 2 oa! s ,io erg s _1 (where P4 = disk pressure 
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Figure 1. ISM gas distribution. Panels a, b, and c display the 
initial gas distribution of the 10 9 , 5x 10 7 , and 10 10 Mq for models 
A, B and C, respectively (see Table 1). 



Figure 2. Total and halo initial column density (cm 2 ) for model 
A. 

in units of (10 4 fc) cm -3 K, ffioo is the disk scale height in 
units of 100 pc and a aj io is the disk sound speed in units 
of 10 kms" 1 ) the superbubble will blowout; thus massive 
starbursts are expected to lead to the superbubble blowout 
phenomena even in massive galaxies such as the Milky Way. 
One can then predict the venting of the hot superbubble 
interior gas through the Rayleigh- Taylor fragmented shell, 
into the extended H I halo where it would push once again 
the outer shock allowing it to build a new shell of swept- 
up halo matter. A coeval starburst with a total mass larger 
than several 10 5 10 6 Mq , also produces an overwhelm- 
ing ionizing photon flux (Fuv > 10 52 photons s" 1 ) and then 
blowout must also allow the leakage of at least a fraction of 
the UV photons emitted by the starburst. These photons are 
here shown to establish an ionized conical Hll region with 
its apex at the starburst. 

It is stellar evolution and the physics of Hll regions 
what defines the time-scales relevant to the evolution, detec- 
tion and impact that massive bursts of star formation may 
have on the ISM. Once the mass of a coeval starburst and 
its IMF are fixed, the integrated mechanical energy power of 
winds and supernovae will remain almost constant for up to 
40 Myr; the life-time of an 8 Mq star. On the other hand, 
the UV photon production rate will remain constant until 
the most massive members of the cluster begin to move away 
from the main sequence (3-4 Myr) to shortly end-up as su- 
pernovae. From then onwards the starburst flux of UV pho- 
tons decays rapidly (as t -5 ; see Beltrametti et al 1982) and 
thus after lOMyr, once the 20 Mq stars have also exploded 
as supernovae, the Hll region gas begins to recombine. De- 
tectability of the ionized gas is thus restricted in the optical 
and radio recombination lines to 10 Myr, while the impact 
of winds and supernovae is a much longer lasting event. One 
can trace in X-rays the thermalized ejected matter either 
within a superbubble or in a galactic wind emanating from 
a nuclear starburst, until it cools either adiabatically or by 
radiation (t coo i — 3fcT/(2nA) ~ 10 8 yr; where n and T are 
the average density and temperature in the bubble, k is the 
Boltzmann constant and A the corresponding cooling rate). 
The remnants of the accelerated ISM can also be detected 
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during a time scale (~ 10 s yr) that well exceeds the energy 
deposition time until they slow down to speeds comparable 
to the random speed of motions of the ISM in their host 
galaxies (~ 10 kms" 1 ). Shells, loops and funnels in a va- 
riety of sizes and shapes have been detected in H I and in 
ionized transitions in many starbursts, dwarf, and irregular 
galaxies (e.g. Marlowe et al. 1995, Hunter & Gallagher 1997 
and Martin 1998) 

Our models are restricted to realistic mechanical energy 
input rates and to the corresponding number of UV photons 
emitted by the massive members in the stellar clusters. Our 
standard model A assumes a 10 6 Mq instantaneous star- 
burst with Z = 0.25 Zq metallicity, and a Salpeter IMF with 
a range of stellar masses between 1 and IOOMq , causing a 
constant mechanical luminosity of 3.5 x 10 40 erg s _1 for up to 
40 Myr. This is consistent with both the synthetic properties 
of starbursts (see Mas-Hesse & Kunth 1991 and Leitherer 
& Heckman 1995) and with the values derived by Marlowe 
et al. (1995) from observations of intermediate mass galax- 
ies with a nuclear starburst. The corresponding number of 
UV photons produced by the massive stars in the cluster 
(N s b = 10 53 photons s _1 ) was regulated to remain constant 
for up to 4Myr, and then, following stellar evolution mod- 
els, to drop as (t/iMyr)~ 5 . A fraction of the UV photons 
produced by the starburst were assumed to be used within 
a starburst region, and therefore only Nuv = £N s b could 
escape out of the central region and ionizes the superbubble 
outer shell, and the extended unperturbed halo. Values of e 
between 0.1 and 0.7 were used in our grid of models. The 
mechanical energy input rate and the UV flux produced by 
the starburst were then linearly scaled with the mass of the 
host galaxy. This led to values of 1.4 xlO 39 erg s _1 and 4 
xlO 51 UV photons s _1 , for the smallest systems, and 3.5 
x 10 42 erg s" 1 and N sb = 10 55 UV photons s" 1 for the giant 
nuclear starbursts in the massive disk galaxies here consid- 
ered. A second calculation of the massive galaxy powered by 
a constant star formation rate instead of a coeval starburst 
is also shown (case C2). In this case the total mechanical en- 
ergy power and UV radiation flux of case CI were assumed 
to remain constant and spread over the 10 8 yr of the con- 
sidered evolution. The main input parameters that define 
our galaxies - starburst strength and boundary conditions - 
are collected in Table 1 which includes in columns 1-3 the 
total mass, the total ISM mass and the central densest ISM 
mass. Columns 4 and 5 list the radius of the dark matter 
distribution and the radius of the galaxy's ISM (see Paper 
I), Chaio and Ccore are the assumed values for the velocity 
dispersion of the dense core and the extended halo compo- 
nents; no is the central number density. Columns 9 and 10 
list the mechanical power of the assumed starburst and its 
total UV output. Column 11 lists the assumed metallicity 
of the ISM and column 12 the starburst life-time. Details 
of the computational approach used to derive the impact of 
the ionizing radiation on the evolving giant remnants and 
the extended halos are given in section 3. 



3 THE NUMERICAL MODEL 

Our hydrodynamical models describe in detail (see Pa- 
per I) the shell morphology, surface density and expan- 
sion velocity, as a function of time. The shock velocity, 



and post-shock densities and temperatures then follow from 
the Rankine-Hugoniot jump conditions. The Mach number 
M s — D a /Cism depends on the ionization state of the ISM 
through the changes in the ambient gas sound velocity 



Cism — 



/ lkT ( 



cff 



1/2 



(3) 



where the mean mass per particle is u n = 14/ limn in the 
neutral region, and fli = 14/23m_ff within the ionized zone. 
The effective gas temperature was taken to be 



Pth 

kn a 



(4) 



where P t h is the thermal pressure and k the Boltzmann's 
constant. We then calculate the post-shock gas adiabatic - 
radiative transition with a tabulated cooling function (Gaetz 
and Salpeter, 1983). 

After blowout, and the disruption of the shell via 
Rayleigh- Taylor instabilities, the hot gas streams out of the 
cavity with a characteristic velocity of the order of it's own 
speed of sound. Therefore we allow the shell to lose mass 
up to the end of acceleration, when the expansion veloc- 
ity equals the sound velocity at the bubble center. After 
blowout the halo gas begins to be collected by a rapid (adi- 
abatic) shock, and a new thin radiative shell forms when the 
characteristic cooling time exceeds the time interval between 
blowout and the current evolution time t. 



3.1 Shell structure 

The structure of the shell of swept-up matter was assumed 
to be able to develop two concentric zones, and the inner 
shell surface was derived from the numerically calculated 
positions. The full thickness of the shell at every Lagrangian 
mesh is 



Is = 



71 s fin 



(5) 



where a is the shell surface density, and fj, n and n s are the 
mean mass per particle and post-shock gas density, respec- 
tively. We then calculate the number of photons required 
to ionize an equivalent spherical shell with the same La- 
grangian segment radius R, surface density, and shell num- 
ber density: 



N c = Y^ R + ls 



r>3i 2 

R \n s a B , 



(6) 



where as = 2.59 x 10 cm s _1 is the H recombination 
coefficient to all levels but the ground state (Osterbrock, 
1989). Taking into account the angle at between the vector 
normal to the shell surface and the radius-vector Ri of a 
particular gridpoint, one can define the ratio 

Ns_ 

N e ' 



P 



(7) 



where N s — Nuv cos on. This allows one to distinguish be- 
tween three obvious cases along each radial direction: a fully 
ionized shell (/? > 1), a neutral shell ((3 = 0), and a partially 
ionized shell (0 < ft < 1), where all UV photons available 
are trapped within the shell. In the later case the condition 



N s = ^-[(R + h) 3 



R 3 ]rvtaB, 



(8) 
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Table 1. Model parameters. 
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1.0 
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0.3 


10 


A2 
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1.4 xlO 3 
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25 
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0.3 


40 


CI 


10.0 
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500.0 


20 


39.6 
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20 


17.0 


350 


10 3 


0.3 


10 


C2 


10.0 


10.0 


500.0 


20 


39.6 


150 


20 


17.0 


140 


4 xlO 2 


0.3 


100 



where U is the thickness of the ionized shell layer should be 
fulfilled. Combining Q with equations ^ and (Q) gives the 
ionized skin thickness: 

U= [p[(R + l s ) s -R A ] + R*] 1/3 ~R. (9) 

The width of the neutral skin then follows from the differ- 
ence between the full shell thickness (||) and the thickness 
of the ionized layer (H) 



l n = R+l s - [P[(R + Is) 3 - R a ] + R 3 ] 



3.2 The ionization of the ISM 



31 1/3 



(10) 



Our simplified approach considers the ionization along a lim- 
ited number of lines of sight (typically along every 15° from 
the axis of symmetry) discarding the shell Doppler shift for 
photon absorption and the UV photons spectral distribu- 
tion. Photons moving along lines of sight that inpinch on 
neutral sections of the shell, are then used up in its ion- 
ization. If the shell segment is in its adiabatic stage, it is 
then transparent to the ionizing radiation for as long as the 
post-shock temperature exceeds 10 4 K. 

The UV photons that escape the shell cause then the 
ionization of the halo. Following an ionization equilibrium 
equation and the on-the-spot approximation, without dust 
absorption (Osterbrock, 1989), the ionization of each shell 
section and the corresponding external interstellar gas is 
considered separately: 

A(R) = N uv - 

( 2 2 f R 2 2 \ 

47r I / n (r)asr dr + n (r)asr dr ) , (11) 
\J Hi JR 2 J 

where Ri and R2 are the inner and outer shell radii respec- 
tively. The shell densities at every cross-point are linearly in- 
terpolated from the values at the nearest Lagrangian points. 
It was also assumed that within the shell the density varies 
linearly a long lines of sight. Then the first integral (Ji) in 
equation (Jill) can be calculated analytically, and reads 



h 



Rl 



nj + (712 



n2^1 

m) 



2m(ri2 — ni) 



Ri 



+ 



(n 2 - ni) 



I N Ri 

m - {ri2 - ni >~f 



{Ri - Ri) 



+ 



(n 2 -m) 2 {Rl-Rf} 



(12) 



where / = R2 — Ri is the shell thickness along a line of sight. 
The second integral in equation ([y]) is calculated numeri- 
cally. 



Stroragren radius 




Figure 3. The trapping of the ionization front. The results of a 
numerical calculation for an ISM gas number density of 10 par- 
ticle cm— 3, and an ionizing photon flux of 10 50 s- 1 . The arrows 
indicate the analytic predictions for the Stromgren radius and the 
characteristic time when all UV photons are trapped within the 
shell of swept-up matter. 



The value of A is first calculated at the galaxy boundary 
(R = Rg)- Positive A values represent the escape of UV 
photons from the galaxy. On the other hand, if A is negative, 
the number of photons is not sufficient to ionize all the halo, 
and the radius of the H II region is then determined by an 
iteration procedure to lie between the shock front and the 
galaxy outer boundary, wherever the value of A becomes 
equal to zero. 

After 3.5 - 4 Myr, the number of UV photons produced 
by the starburst rapidly drops. Also, between 4-13 Myr 
(Models A and B) the leading shock wave moving in the 
halo becomes radiative and the density within the shell dras- 
tically increases (as the square of the Mach number) , causing 
a larger recombination rate. Both proceses reduce the UV 
flux available to keep the galaxy halo fully ionized, and lead 
to the trapping of the ionization front within the expanding 
shell. 

Our procedure to estimate the trapping of the ioniza- 
tion front, was checked against the analytic formula for a 
spherical bubble expanding in the preionized homogeneous 
medium (Comeron 1997), and we found our numerical re- 
sults in full agreement with analytic predictions (see Fig. ^). 

Our procedure also considers the fact that upon pho- 
toionization, the recombination time (r rec = l/n e QA ~ 
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10 5 /rihaio yr) of the low density halo gas may exceed 
the characteristic dynamical time of the superbubble. This 
causes the ionized halo to become temporarily transparent 
to the ionizing radiation. We thus calculate the character- 
istic recombination time r rec along different lines of sight 
and compare this value with a time interval r = t — tbiowout 
between a current time t and the time when a shell section 
becomes transparent to the starburst UV radiation, causing 
the ionization of the halo. The number of ionizing photons 
which escape the galaxy in any direction follows from the 
equation 

N esc (Q) = N UV - Ntrap(Rr), (13) 

where G is the angle between the direction considered and 
the symmetry axis. If T rec (Rc) < T, the halo can only be 
kept ionized if there is a continuous input of UV photons, 
and the radius R T equals the halo radius (R T = Kg). Oth- 
erwise, if Trec(RG) > t, recombination within the densest 
sections of the galaxy halo will diminish the escape of UV 
photons. In this case the radius R T follows from the condi- 
tion T rec (R T ) = T. 



4 A TYPICAL EVOLUTIONARY SEQUENCE 

Fig. shows the evolutionary sequence of our basic model 
Al. It compares the remnant expansion speed along different 
angles (O) from the symmetry axis with the galaxy escape 
speed, and indicates the maximum size acquired by the su- 
perbubble as a function of time. The shell expansion velocity 
shows at first a continuous deceleration, followed by a rapid 
acceleration of the pole sections of the remnant immediately 
after blowout and shell disruption. During the acceleration 
phase the shock builds a new shell of swept up halo matter. 
Afterwards, the remnant speed begins to decline even along 
the direction perpendicular to the galaxy plane (O = 0), as 
the shock gathers more halo matter. In this way, the rem- 
nant ends up expanding with speeds well below the galaxy 
escape velocity V esc (see also paper I). Nevertheless, after 40 
Myrs (the starburst lifetime) a huge remnant has developed 
and extends in the pole direction almost 5 kpc out of the 
starburst region (see Fig. ^|a). This late phase in the evolu- 
tion of the remnant, as stated earlier, is unobservable in the 
optical regime. 

Fig. p»a-f) show the growth of the superbubble as well 
as the development of the extended ionized zone and it's 
evolution as a function of time. Before 1.5 Myr the dense 
shell traps all UV photons escaping the central starburst re- 
gion. However, after blowout and shell disruption, the UV 
photons escaping the central H II region produce an initially 
narrow cone of ionized interstellar gas around the symme- 
try axis. The conical Hll region becomes rapidly broader, 
and at t ~ 2.5 Myr it reaches it's maximum opening angle. 
Note that the low density halo matter has a long recom- 
bination time {t rec = (an) -1 ~ 10 5 yr/n; where a is the 
recombination coefficient and n the local density) and thus 
upon photoionization it becomes undetectable both in the 
optical regime and in the 21 cm line and it also becomes 
transparent to further UV radiation arising from the star- 
burst. After 5 Myr the superbubble leading shock makes a 
transition from the adiabatic to the radiative regime and 
matter in the swept up shell begins to recombine. Given the 



reduced number of ionizing photons produced by the star- 
burst after 4 Myr and the dilution factor experienced by 
the UV radiation to reach the recombining shell, rapidly all 
available photons begin to be used within the shell inhibit- 
ing the further leakage of UV radiation into the halo and 
out of the galaxy. Then, after 8 Myr, the densest parts of 
the still unshocked halo begin to recombine. Recombination 
proceeds outwards from the densest unshocked halo regions, 
forming an increasingly wider layer of neutral gas ahead of 
the expanding shell (see Fig. ^c). 

It is not a straightforward issue to estimate the frac- 
tion of Lyman continuum photons that leaks out of galaxies. 
The only direct observations of the Lyman break in nearby 
star-forming galaxies are those made by the Hopkins UV 
telescope (HUT) (Leitherer et al. 1995). The analysis sug- 
gests that less than 3 percent of the intrinsic Lyman con- 
tinuum photons escape into the intergalactic medium (Lei- 
therer et al. 1995). A more recent analysis of the same data 
indicates less restrictive values, between 3 percent and 57 
percent (Hurwitz et al. 1997). Our models predict an early 
evolutionary phase (between blowout and the trapping of 
the ionization front within the radiative expanding shell i. e. 
between 3 Myr and 5 Myr for our Al model) during which a 
large amount of the UV radiation could leak out of a galaxy 
into the intergalactic medium (see Fig. |i]a). Our predictions 
seem to be well in line with the observations given that all 
four star-bursts observed with HUT are much older than 
5 Myr (Gonzalez-Delgado et al. 1998) i.e. a time when no 
escape is expected. 

The models have assumed that the production of Ly- 
man continuum radiation remains constant for up to 4 xlO 6 
years and then, as the ionizing cluster ages, it drops as t -5 
causing after 10 7 years (the expected life-time of the central 
H II region) a depletion in the ionizing radiation of more than 
two orders of magnitude. An important fraction (e = 0.7) of 
the initial UV photon production rate was assumed, in all 
cases shown in Fig. [| to leak out of the central Hll region 
causing the ionization of the halo and of the recombining 
sections of the expanding shell. Fig. ^ shows the fraction 
of UV photons able to escape the galaxy along the coni- 
cal H II region, as well as the time interval during which this 
becomes possible. Our results, the same as those from calcu- 
lations with other values of e, with the exception of the case 
that assumes a constant star formation rate (Fig. ^d), all 
show that photoionization of the halo and the escape of UV 
radiation from the galaxy occurs regardless of the produc- 
tion of UV radiation (and the assumed value of e) and in all 
cases limited to 5 x 10 6 years. Given the importance of the 
escape of the Lyman continuum radiation in relation with 
estimates of star formation rates, a more detailed analysis 
of this point is in preparation. 

The rapid changes in ionization of the galaxies' ISM 
have an impact on the neutral gas column density arising 
from the expanding shell and the halo, and these are central 
to the transport of the Lya line emitted by the nuclear H II 
region which otherwise will be completely absorbed within 
the extended Hi halo. If one for example considers the di- 
rection along the axis of symmetry, one would see that im- 
mediately after blowout, the ionization of the low density 
halo, not only favours the direct escape of UV photons out 
of the galaxy (see Fig. ^i) but also generates a transparent 



© 1994 RAS, MNRAS 000, 



Superbubbles in star-forming galaxies 7 



BOO 




' ' ' 1 ' 


,,,,,,,, 


- 


700 


: 


= 0° 

/ 






600 


: 








T m 500 










km 










400 
















Slop 




i 300 




1 ; 9 = 30° 
\V ( = 60° 






100 





\ / 










= 90° 

, , I , , , , 1 





20 

Time (Myr) 




10 20 30 40 50 60 70 
Time (Myr) 





10 20 30 

Time (Myr) 



10 20 30 

Time (Myr) 



Figure 4. The evolution of superbubbles. Panels a-d display the results for models Al, B, CI and C2. Solid lines represent the expansion 
velocities of the superbubbles along different lines of sight (0), and as a function of time. These are to be compared with the dotted line 
which indicates the local escape velocity from a galaxy center at a distance comparable to the maximum dimension of the superbubble 
(along © = 0), indicated by the dashed lines. 



region for the Lya photons produced upon recombination in 
the central Hll region (see Fig. [j]). 

On the other hand, the leading shock transition from 
adiabatic to radiative is to deplete the number of UV pho- 
tons into the halo and out of the galaxy and eventually to 
cause the trapping of the ionization front within the ex- 
panding shell. This situation favours the development of an 
expanding H I layer and leads to the appropriate conditions 
(column density and velocity) for the appearance of a blue 
shifted Lya absorption feature. In model Al, full ionization 
of the halo is completed before 3 Myr, and from then on- 
wards and until the expanding shell becomes radiative and 
traps the ionization front, an observer looking at the central 
Hll region within the ionized halo cone should see the full 
Lya line in emission. The trapping of the ionization front 
occurs at t ~ 5 Myrs, when the central starburst is still pro- 
ducing the sufficient number of UV photons to maintain the 
ionization of the central H II region. This indicates a high 
probability to observe a blue shifted Lya absorption feature 
in intermediate mass galaxies. Recombination in the halo 
starts after 9 Myr, rapidly enhancing the column density of 
neutral material and with it the full absorption of the Lya 
line produced in the central Hll region, irrespective of the 
line of sight into the galaxy. Panels a, b and c in Fig. R show 



the run of the H I column density along the axis of symme- 
try for the standard model A with three different metalici- 
ties (models A1-A3). The change in the assumed ISM metal 
content simply speeds (or delays) the onset of radiative cool- 
ing in the expanding shell, and with it the trapping of the 
ionization front within the shell. Consequently, superbub- 
bles evolving in more metal rich galaxies (as in model A2) 
will have a shorter lasting phase in which to show Lya in 
emission, and a correspondingly longer phase causing a blue- 
shifted Lya absorption feature. The converse also holds for 
more metal poor galaxies (as in model A3). 

Despite the similar evolution of the shell expansion ve- 
locity, which in all cases shows the same behaviour (at first 
a continuous deceleration, followed by a rapid acceleration 
phase after blowout, and a strong deceleration as the bub- 
ble ploughs through the extended H I halo) the evolution of 
superbubbles in other galaxies departs in a variety of ways 
from our standard case. The fact that remnants of powerful 
starbursts in massive galaxies, with steep density gradients, 
retain high expansion speeds (of several hundreds of km s _1 ) 
along the polar sectors of their superbubbles for longer pe- 
riods than their counterparts in less massive galaxies (see 
Fig. ^]) is to be kept in mind. The temperature acquired by 
the shocked gas being directly related to the shock speed 
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Figure 5. The extended conical HII regions. Cross-sectional views along the symmetry axis display the growth of the superbubblc (solid 
contours) and the development of the conical HII region (shaded area) in model Al (panels a-c). Note that the HII region extends to 
the edge of the galaxy, reducing in this way the column density of neutral material within a wide cone. Panels a-c display the evolution 
after 3.1, 3.5 and 15 Myr. Panel c shows the growth of an increasingly wider shell of neutral matter generated by recombinations in the 
halo at late evolutionary times, d-f show the results after 15 Myr of evolution for case B, CI, and C2. Note that the low densities in the 
halos of these galaxies have inhibited recombinations even after the main H II region life-time. 



{Tshock = 1.4 x 10 7 (n s h oc fe/1000kms -1 ) 2 ) implies that these 
shells are to become radiative at a time that well exceeds 
both the Hll region life-time (~ lOMyr) and the 5 xlO 7 
years span during which the starburst deposits its mechani- 
cal energy. Therefore, these galaxies are not likely to produce 
a blue shiffted Lya absorption feature. 

Low mass galaxies appear at first glance as the best 
candidates to lose their starburst processed metals, or even 
their ISM (Fig. [jjb). However, these systems also experience 
a rapid shock deceleration as their superbubbles evolve into 
the galaxy halos, leading as in the standard case A to speeds 
well below the galaxy escape velocity (see Fig. ^>). Also, as 
shown in Paper I, if the mass of the galaxy is initially further 
centrally condensed, then the larger densities and stronger 
cooling imply an amount of energy much larger than what 
is observed (see Marlowe et al. 1995) for the superbubble to 
reach the edge of the galaxy. On the other hand, massive 
disk galaxies with a powerful nuclear starburst generate gi- 
ant superbubbles able to reach supersonically the outskirts 
of the galaxy, dumping their processed matter into the inter- 
galactic medium. This happens both in the coeval starburst 
model (CI) and in the constant star formation rate one (C2) 
at a time that well surpasses 40Myr. 



Fig. |H|i-f show a cross-sectional view of the superbubble 
and the conical Hll region developed after 15Myr of evolu- 
tion for cases Bl, CI and C2, respectively. The low density 
halos in these galaxies inhibit recombination for times that 
well surpase the H II region life-time. These panels are to be 
compared with Fig. ^c where the recombination of the halo 
matter (after 15Myr of evolution) has led to a broad shell 
of recombined matter ahead of the evolving superbubble in 
an intermediate mass galaxy (case Al). 

Strong radiative cooling within the expanding shell pro- 
motes its transition from adiabatic to radiative and this - if 
it occurs before 5 x 10 6 yr (see Fig. []) - has an important im- 
pact on the number of UV photons leaking into the halo and 
out of the galaxy. In both, the low and high mass systems 
here considered, the radiative shell forms much later than in 
the intermediate mass galaxies, long after the massive stars 
have moved off the main sequence. 

Despite the large value of e = 0.7 used in all calculations 
presented here, in all models with a coeval starburst the 
leakage of UV photons into the halo and possibly out of the 
galaxy is restricted to 5Myr of evolution by the boundary 
condition imposed by stellar evolution (Fuv ~ ^~ 5 )- The 
drastic drop in the photon production rate after the most 
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Figure 6. The number of UV photons that escape from a galaxy in models A, B, CI and C2 are shown in panels a,b,c and d respectively, 
as a function of time. The escape of UV radiation does not start until the whole galaxy halo is ionized, and in all coeval starburst models 
(A— CI) it is restricted to 5Myr by the condition imposed by stellar evolution. Note that given the low densities in the galaxy halos and 
the corresponding recombination times, these remain ionized for times that well exceed the HII region life-time (10 Myr). 



massive stars begin to move off the main sequence, inhibits 
UV photons to even leak out of the central H II egion, and 
thus clearly calculations with smaller values of e display the 
same upper limit (see Fig. ^|a-c). The situation is different in 
the constant star formation rate (case C2) where the number 
of UV photons escaping the galaxy remains almost constant 
throughout the evolution. This is however also a result of 
the long recombination time in the low density halo which 
makes it transparent to the starburst UV flux. 



5 DISCUSSION 

Our numerical calculations of superbubbles evolving in gas- 
rich galaxies, accounting for the mechanical energy deposi- 
tion as well as for the ionization produced by the massive 
members of a starburst, have led to a more complete pic- 
ture of the impact that massive star formation may have on 
a galaxy ISM. We have centered our attention on the time 
evolution of superbubbles to evaluate the possible burst of 
processed matter into the intergalactic medium and found 
that this is only likely for powerful nuclear bursts in disk 
galaxies. Small and intermediate mass galaxies host super- 
bubbles unable to reach their galaxy outskirts. Furthermore, 
one should note that the first lOMyr of evolution, i.e. the 



H II region life-time, is only a small fraction of the time re- 
quired for superbubbles to grow and reach their largest di- 
mensions, and thus studies in optical wavelengths can only 
address the beginning of the evolution and the possible burst 
into the intergalactic space occurs at much latter times. We 
have also looked at the extended conical Hll region gener- 
ated by the escape of UV photons from the central starburst 
once blowout of the superbubble into the extended H I halo 
takes place. The early blowout into the extended H I halo is 
promoted by the steep density drop away from the galaxy 
plane, causing the burst of the thermalized ejected matter 
as well as the leak of UV radiation into the extended halo 
(see Fig. |^a, b). The evolution of the Hll region is central 
to the transport and detection of Lycv emission from the 
central Hll region. For example, once recombination starts 
in the fast expanding shell, it will cause a correspondingly 
blueshifted Lya emission, as depicted in Fig. |^c. Once the 
shell presents a large column density (~ 10 19 cm -2 ), as it 
grows to dimensions of a few kpc, it will trap the ioniza- 
tion front. This is promoted by the large shell densities and 
the geometrical dilution of the ionizing radiation. Note that 
from then onwards, recombinations in the shell will inhibit 
the further escape of ionizing photons from the galaxy (com- 
pare Fig. 0b, c, and d). The trapping of the ionization front, 
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Figure 7. The neutral gas column density variation along the 
symmetry axis for models A1,A2, A3, and B are shown respec- 
tively in panels a — d. Solid lines indicate the unshocked com- 
ponent and the dotted lines the column density arising from the 
recombining shell. 



makes the shell acquire a multiple structure with a photoion- 
ized inner edge, a steadily growing central zone of H I, and an 
outer collisionally ionized sector where the recently shocked 
ionized halo gas is steadily incorporated. The growth of the 
central layer eventually will cause the sufficient scattering 
and absorption of the Lya photons emitted by the central 
Hll region, leading to a blueshifted Lya absorption. Note 
also that for as long as recombinations continue to occur at 
the leading edge of the shell, a blueshifted Lya in emission 
will appear superposed on the blueshifted absorption feature 
(see Fig. |^d). Recombination at the leading edge of the shell 
will become steadily less frequent, depleting the blueshifted 
Lya in emission. This is due when the shell and its lead- 
ing shock move into the outer less dense regions of the halo, 
and the shell recombination time, despite the compression at 
the shock, becomes larger than the dynamical time. At this 
stage, an observer looking along the conical Hll region will 
detect a P-Cygni-like Lya line profile as shown in Fig. ^e. 

The geometrical dilution of the UV flux will begin to 
make an impact as the superbubble grows large. This, and 
the drop in the UV photon production rate caused by the 
death of the most massive stars after t — t ms , will enhance 
the column density of neutral material in the central zone of 
the recombined shell to eventually cause the full saturated 
absorption of the Lya line (see Fig. ^f). Full saturated ab- 
sorption has usually been accounted for by the large column 



density of the extended Hi envelope of these galaxies and 
thus, as in all models, many different orientations will suit 
the observations. In our scenario however, also when observ- 
ing within the solid angle defined by the conical H II region 
formed after blowout, such broad absorptions could arise 
well after blowout, once a large column density of shocked 
neutral material (N > 10 20 cm -2 ) has formed ahead of the 
trapped ionization front. 

The solid angle subtended by the ionized cone will be 
a rapidly changing function of time, particularly during the 
early stages of evolution, immediately after blowout. How- 
ever, numerical experiments and observations (see Paper I, 
and Tenorio-Tagle & Munoz-Tunon 1997, 1998 and refer- 
ences therein) restrict this to a maximum value of about 
70° , with the walls of the superbubble near the galaxy plane 
inhibiting its further growth. 

P-Cygni Lya profiles are predicted when observing 
along the angle subtended by the conical Hll region but 
only once the ionization front is trapped by the sector of 
the superbubble shell evolving into the extended halo. This 
will produce the fast moving layer of Hi at the superbub- 
ble shell, here thought to be responsible for the partial ab- 
sorption observed in sources such as Haro 2, ESO 400-G043 
(which probably exhibits a secondary blueshifted Lya emis- 
sion) and ESO 350-IG038 (Kunth et al. 1998). In the latter 
case however, the profile is not typical of a clear P-Cygni 
profile. Instead the underlying damped Lya absorption ex- 
tends beyond the red of the line emission. Damped Lya ab- 
sorption is shown by several galaxies. We note that these 
objects are all gas rich dwarf galaxies whereas in most cases 
but Haro 2, the galaxies that exhibit Lya in emission or with 
a P-Cygni profile, are on the higher luminosity side of the 
distribution (M < —18). Pure Lya emission is observed in 
C0840+1201 and T1247-232 (Terlevich et al. 1993; IUE) or 
T1214-277 (Thuan & Izotov 1997; HST). Such a line implies 
no absorption and thus no H I gas between the starburst H II 
region and the observer, as when observing the central Hll 
region after the superbubble blowout, within the conical Hll 
region sector carved in the extended Hi halo. 

The main implication of our calculations is that it is the 
feedback from the massive stars what - through ionization 
and the evolution of superbubbles - leads to the large variety 
of Lya emission profiles. The escape of Lya photons depends 
sensitively on the column density of the neutral gas and dust 
following the suggestion that the attenuation by dust is en- 
hanced by scattering with hydrogen atoms. Note that apart 
from the observed profiles: Lya in emission, P-Cygni-like 
profiles and full saturated absorption (as in Fig. H a, b, e 
and f) the scenario predicts also secondary blueshifted Lya 
emission profiles emanatingfrom the rapidly expanding and 
recombining shell (see Fig. pp and d). If massive star forma- 
tion leads also to networks of shells such as those observed 
in 30 Dor (Chu & Kennicutt 1992) one should also expect a 
forest of Lya in emission arising from recombinations in the 
various expanding shells in the network. 
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Figure 8. The evolution of superbubbles and the hya emission profile expected from star-forming galaxies. The panels present the 
sequence of changes in the structure of the ISM produced by the evolution of a massive starburst. On the right hand-side of each panel, 
the hya profile expected is shown and the scale corresponds to model Al at times (a)< 1.5Myr, (b)~2.5Myr, (c)~4Myr, (d)~5Myr, 
(e)~6Myr and (f)>8Myr. The sequence starts with a giant superbubble which is about to burst into the extended Hi halo of a galaxy, 
as its shell of swept up matter becomes Rayleigh - Taylor unstable (a). The instability leads to the rupture of the shell and with it to 
the exit of the thermalized stellar ejecta which then freely streams and promotes the development of a new shock and a new shell of 
swept up halo material (b). At the same time, the stellar uv radiation is able to leak out into the halo, and even beyond the galaxy outer 
edge, leaving behind an extended, low density, conical H II region. Recombinations in the shell of swept up halo matter should become 
increasingly more important as a function of time, causing a blueshifted hya emission and a depletion in the number of ionizing photons 
escaping the galaxy. The latter will be totally depleted once the ionization front is trapped by the recombining shell (c and d). The last 
stages in the evolution enhance the column density of H I in the recombining shell leading to a P-Cygni hya profile and eventually to a 
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